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ABSTRACT 

We explore possible cosmological consequences of a running Newton's constant as suggested 

by the non-trivial ultraviolet fixed point scenario in the quantum field-theoretic treatment of Ein- 
stein gravity with a cosmological constant term. In particular we focus here on what possible 
effects the scale-dependent coupling might have on large scale cosmological density perturbations. 
Starting from a set of manifestly covariant effective field equations derived earlier, we systemat- 
ically develop the linear theory of density perturbations for a non-relativistic, pressure-less fluid. 
The result is a modified equation for the matter density contrast, which can be solved and thus 
provides an estimate for the growth index parameter 7 in the presence of a running G. We complete 
our analysis by comparing the fully relativistic treatment with the corresponding results for the 
non-relativistic (Newtonian) case, the latter also with a weakly scale dependent G. 
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1 Introduction 



Recent years have seen the development of a bewildering variety of alternative theories of gravity, 
in addition to the more traditional alternate theories, which used to include Brans-Dicke, tensor- 
scalar, tensor-vector-scalar, higher derivative, effective quantum gravity and supergravity theories. 
Some of the new additions to the already rather long list include dilaton gravity, f{R) and f{G) 
gravity, Chern-Simons gravity, conformal gravity, torsion gravity, loop quantum gravity, holographic 
modified gravity, MOG gravity, asymmetric brane gravity, massive gravity and minimally modified 
self-dual gravity, just to cite a few representative examples. All of these theories eventually predict 
some level of deviation from classical gravity, which is often parametrized either by a suitable set 
post-Newtonian parameters, or more recently by the introduction of a slip function [1, 2]. The 
latter has been quite useful in describing deviations from classical GR, and specifically from the 
standard kCMD model, when analyzing the latest cosmological CMB, weak lensing, supernovae 
and galaxy clustering data. 

In this paper we will focus on the analysis of departures from GR in the growth history of 
matter perturbations, within the narrow context of the non-trivial ultraviolet fixed point scenario 
for Einstein gravity with a cosmological term. Thus instead of looking at deviations from GR 
at very short distances, due to new interactions such as the ones suggested by string theories 
[3], we will be considering here infrared effects, which could therefore become manifest at very 
large distances. The classical theory of small density perturbations is by now well developed in 
standard textbooks, and the resulting theoretical predictions for the growth exponents are simple 
to state, and well understood. Except possibly on the very largest scales, where the data so far is 
still rather limited, the predictions agree quite well with current astrophysical observations. Here 
we will be interested in computing and predicting possible small deviations in the growth history 
of matter perturbations, and specifically in the values of the growth exponents, arising from a 
very specific scenario, namely a weakly scale-dependent gravitational coupling, whose value very 
gradually increases with distance. 

The specific nature of the scenario we will be investigating here is motivated by the treatment 
of field-theoretic models of quantum gravity, based on the Einstein action with a bare cosmological 
term. Its long distance scaling properties are derived from the existence of a non-trivial ultraviolet 
fixed point of the renormalization group in Newton's constant G. The latter is inaccessible by direct 
perturbation theory in four dimensions, and can be shown to radically alter the short- and long- 
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distance behavior of the theory when compared to more naive expectations. The renormalization 
group origin of such fixed points was first discussed in detail by Wilson for scalar and self-coupled 
fermion theories [4] . The general field theoretic methods were later extended and applied to gravity, 
where they are now referred to as the non-trivial fixed point scenario or asymptotic safety [5] . It is 
fair to say that so far this is the only field-theoretic approach known to work consistently in other 
not perturbatively renormalizable theories, such as the non-linear sigma model. While perhaps 
still a bit mundane in the context of gravity, such non-trivial fixed points are well studied and 
well understood in statistical field theory, where they generally describe phase transitions between 
ordered and disordered ground states, or between weakly coupled and condensed states. 

The paper is organized as follows. First we recall the effective covariant field equations describing 
the running of G, and describe the nature of the objects and parameters entering the quantum 
non-local corrections. We then discuss the zeroth order (in the fluctuations) field equations and 
energy-momentum conservation equations for the standard homogeneous isotropic metric, with a 
running G. Later we extend the formalism to deal with small metric and matter perturbations, 
and derive the relevant field and energy conservation equations to first order in the perturbations. 
After showing the overall consistency of the derived equations, we proceed to derive the modified 
differential equation for the density contrast 6{t). Later this is re- written, following customary 
procedures, as a function of the scale factor as 6 {a). The resulting differential equation for the 
density contrast is then solved and the results for the growth exponents compared to the standard 
classical result. The conclusions provide an interpretation of the theoretical results and their 
associated uncertainties vis-a-vis present and future high precision galaxy clustering measurements. 



2 Running Newton's Constant G(n) 



Originally the running of G was computed cither on the lattice directly in four dimensions 
[6, 7, 8], or in the continuum within the framework of the background field expansion applied to 
2 + € spacetime dimensions [5, 9] and later using truncation methods applied in 4d [10]. In either 
case one obtains a momentum dependent G{k'^), which needs to be eventually re-expressed in a 
coordinate-independent way, so that it can be usefully applied to more general problems involving 
arbitrary background geometries. 

The first step in analyzing the consequences of a running of G is therefore to re-write the 
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= Go 



1 + Co 



(2.1) 



expression for G{k'^) in a coordinate-independent way, either by the use of a non-local Vilkovisky- 
type effective gravity action [11, 12], or by the use of a set of consistent effective field equations. 
In going from momentum to position space one usually employs fc^ — □, which then gives for 
the quantum-mechanical running of the gravitational coupling the replacement G — > G(n). One 
then finds that the running of G is given in the vicinity of the UV fixed point by 

where □ = g^'^V ^iVv is the covariant d'Alembertian, and the dots represent higher order terms in 
an expansion in l/(^^n). Current evidence from Euclidean lattice quantum gravity points toward 
Co > (implying infrared growth) and v —\ [8]. 

Within the quantum-field-theoretic renormalization group treatment, this last quantity arises as 
the integration constant of the Callan-Symanzik renormalization group equations. One challenging 
issue therefore, and of great relevance to the physical interpretation of the results, is a correct 
identification of the renormalization group invariant scale ^. A number of arguments can be given 
(see below) in support of the suggestion that the infrared scale ^ (very much analogous to the ^jjg 
of QCD) can in fact be very large, even cosmological, in the gravity case. Prom these arguments 
one would then infer that the constant Go can, to a very close approximation, be identified with 
the laboratory value of Newton's constant, \/Go ~ 1.6 x lO^^^cm. 

The appearance of the d'Alembertian □ in the running of G naturally leads to both a non-local 
effective gravitational action, and a corresponding set of non-local modified field equations. Instead 
of the ordinary Einstein field equations with constant G 

Rixv - k9nuR + Xgt^u = SttGT^^ , (2.2) 

one is now lead to consider the modified effective field equations 

R^u - \gt.uR + Xg^^v = Stt G(n) T^, (2.3) 

with a new non-local term due to the G{0). By being manifestly covariant they still satisfy some 
of the basic requirements for a set of consistent field equations incorporating the running of G. 
Not unexpectedly though, the new nonlocal equations are much harder to solve than the original 
classical field equations for constant G. 

It is instructive to note, as already pointed out in [13], that the effective non-local field equations 
of Eq. (2.3) can be re-cast in a form very similar to the classical field equations, but with a new 
source term T^^i, = [G(n)/Go] T^i^ defined as the effective, or gravitationally dressed, energy- 
momentum tensor. Ultimately the consistency of the effective field equations demands that it be 
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exactly conserved, in consideration of the contracted Bianchi identity satisfied by the Ricci tensor. 
In this picture, therefore, the running of G can be viewed as contributing to a sort of a vacuum 
fluid, introduced in order to account for the new gravitational vacuum polarization contribution. 

More on the technical side, and mainly due the appearance of a negative fractional exponent in 
Eq. (2.1) , the covariant operator appearing in the expression for G(n) has to be suitably defined 
by analytic continuation. This can be done, for example, by computing for positive integer 
n, and then analytically continuing to n — >■ —1/21/ [13]. Equivalently, (?(□) can be defined via a 
suitable regulated parametric integral representation [16], such as 



As far as the calculations in this paper are concerned, it will not be necessary to commit oneself 
to an unduly specific form for the running of G{0). Thus for example, although the lattice gravity 
results only allow for a non-degenerate phase for the case cq > 0, it will nevertheless be possible 
later to have either sign for the correction in Eq. (2.1), in the sense that the very existence of a non- 
trivial ultraviolet fixed point implies in principle the appearance of two physically distinct phases, 
each of which might or might not be physically realized due to issues of non-perturbative stability. 
Observation could then be used, in principle, to constrain one or the other choice. Furthermore, 
the value of the exponent need not to be specified until the very end of the calculation, so that 
most of the results can be kept general. ^ 

The situation regarding the running of G is perhaps most easily illustrated close and above two 
dimensions, where the gravitational coupling becomes dimensionless, G ^ with A the ultra- 
violet cutoff required to regularize the theory (a similar and completely parallel line of arguments 
and results can in fact be presented for the 4d lattice theory as well, but a discussion of renormal- 
ization on the lattice ends up being inevitably quite a bit less transparent [8, 6]). There the theory 
appears perturbatively renormalizable, so that the full machinery of covariant renormalization and 
of the renormalization group can in principle be applied, following Wilson's dimensional expansion 
method, now formulated as a double expansion in G and e = d — 2 [5, 9]. Both here and on the 
lattice a renormalization of the bare cosmological constant, besides being gauge-dependent, is also 
physically meaningless, as it can be reabsorbed by a trivial rescaling of the metric; the latter is 

running cosmological constant A(fc) A(n) causes a number of mathematical inconsistencies [13] within the 
manifestly covariant framework, described here by the effective field equations of Eq. (2.3). Indeed if one assumes 
that A(D) ~ (^^□)~'', where a is a (positive or negative) power, then for example the infrared regulated expression 
in Eq. (2.4) gives no running of A, after using Va<7,i^ = 0. This last conclusion is in agreement with the field-theoretic 
results of the non-trivial renormalization group fixed point scenario, thereby providing perhaps an independent 
consistency check. 




(2.4) 
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needed in order to recover the proper normalization of the vohime term in the path integral, thus 
avoiding spurious renormalization effects, as discussed in [6, 8, 9]. 

In momentum space the result corresponding to Eq. (2.1), and allowing now possibly for either 
sign in front of the correction, is 



1 \l/2!^ 

l±^o(^j +••• 



(2.5) 



with Co a positive constant, and ^ the new, genuinely nonperturbativc, gravity scale. ^ Consequently 
the above expression for G{k^) can be used whenever the full generality of the manifestly covariant 
expression in Eq. (2.1) is not really needed, for example when dealing with the Newtonian (non- 
relativistic) limit. 

The choice of + or — sign is ultimately determined from whether one is initially to the left (+), 
or to right (-) of the fixed point Go, in which case the effective G(fc^) decreases or, respectively, 
increases as one flows away from the ultraviolet fixed point towards lower momenta, or larger 
distances. Physically the two solutions represent of course gravitational screening (G < Go) or 
anti-screening (G > Go). 

It is crucial that the quantum correction involves a new physical, renormalization group invari- 
ant, scale ^, whose value cannot be fixed by a perturbative calculation, and whose absolute size 
determines the comparison scale for the new non-local quantum effects. It should therefore be right- 
fully considered as the gravity analog of the celebrated gauge theory scaling violation parameter 
^Ws- terms of the bare gravitational coupling G(A) it is given by 

where /3(G) is the Callan-Symanzik beta function for G (which can be given explicitly, for example, 
in the 2+e expansion to a given loop order, or can be computed on the lattice). It is then more or less 
a direct consequence of the renormalization group that the value of the constant determines the 
coefficient cq in Eq. (2.1), cq = 1/{A^^'^ Gq). The non-perturbative lattice formulation of quantum 
gravity then allows an explicit and direct computation of ^4^, and therefore of the coefficient cq in 
G(n) [8. 6]. 



*A properly infrared regulated version of the above expression, here with the choice of -|- sign, would read 

-2 \ 1/21' 



G(fe") ~ Go 



1 + ^=0(^^1 +.. 



(2.6) 



Then for large distances r S> C the gravitational coupling no longer exhibits the spurious infrared divergence, but 
instead approaches the finite value Goo ~ (1 -|- co -|- . . .) Go- 
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Physically it would seem at first that the non-perturbative scale ^ could take any value, including 
a very small one - based on the naive estimate ^ ~ /p - which would then of course preclude 
any observable quantum effects in the foreseeable future. But a number of recent results for the 
gravitational Wilson loop on the Euclidean lattice at strong coupling, giving an area law, and 
their subsequent interpretation in light of the observed large scale semiclassical curvature [15], 
would suggest otherwise, namely that the non-perturbative scale ^ appears in fact to be related to 
macroscopic curvature. Prom astrophysical observation the average curvature on very large scales, 
or, stated in somewhat better terms, the measured cosmological constant A, is very small. This 
would then suggest that the new scale ^ can be very large, even cosmological, 

which would then give a more concrete quantitative estimate for the scale in the G(n) of Eq. (2.1), 
namely $, ~ l/^/\J^ ~ 1.51 x 10^^ cm. Indeed for quantum gravity no other suitable infrared cutoff 
presents itself, so that A can almost be considered as the only "natural" candidate to take on the 
role of a (generally covariant) infrared regulator or graviton mass-like parameter. 

Finally let us mention here briefly and for completeness that for a limited number of metrics 
it has been possible, after some considerable work, to find exact solutions, in some regime, to the 
above effective nonlocal field equations. One such case is the static isotropic metric, where in the 
limit r ^ 2AIG one can obtain an explicit solution for the metric coefficients A(r) = 1/B{r), 
leading eventually to the rather simple result [14] 



with m = consistent with a gradual slow increase of G{r) with distance. One amusing aspect 
of the exact solution in the static isotropic case is that no consistent solution can be found unless 
= 1/3 exactly in four dimensions, and similarly v = l/(d — 1) vn dimensions c? > 4 [14], lending 
further support, and independently of the lattice theory results, to this particular value for u in 
four dimensions. 

®We have pointed out before that the result for G(r) is in a number of ways reminiscent of the analogous QED 
result (known as the Uehling correction to the Coulomb potential in atoms) 



In the gravity case the correction is not a log but a power, which is what one would naively expect from a perturbatively 
non-renormalizable theory. In gravity, the infrared cutoff due in QED to the finite physical electron mass is naturally 
replaced by the physical cosmological constant; the magnitude of neither one of these two quantities can be predicted 
by the fundamental theory. 




(2.9) 




(2.10) 
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2.1 (Zeroth Order) Effective Field Equations with G(n) 



A scale dependent Newton's constant is expected to lead to small modifications of the standard 
cosmological solutions to the Einstein field equations. Here we will summarize what modifications 
are expected from the effective field equations on the basis of G(n), as given in Eq. (2.1), which itself 
originates in Eq. (2.5). The starting point are the quantum effective field equations of Eq. (2.3), 
with G{0) defined in Eq. (2.1). In the Fricdmann-Lemaitre-Robertson- Walker (FLRW) framework 
these are applied to the standard homogeneous isotropic metric 

dr^ = dt2-a2(i) 1^^^ +r2 (ci^^^sin^^ A; = 0,±1. (2.11) 

In the following we will mainly consider the case A; = (spatially flat universe). It should be noted 
that there are in fact two related quantum contributions to the effective covariant field equations. 
The first one arises because of the presence of a non- vanishing cosmological constant A 2± 3/^^, 
caused by the non-perturbative vacuum condensate. As in the case of standard FLRW cosmology, 
this is expected to be the dominant contributions at large times t, and gives an exponential (for 
A > 0), or cyclic (for A < 0) expansion of the scale factor. The second contribution arises because of 
the explicit running of G{0) in the effective field equations. The next step therefore is a systematic 
examination of the nature of the solutions to the full effective field equations, with G(n) involving 
the relevant covariant d'Alembertian operator 

□ = 5^^V^V. (2.12) 

acting on second rank tensors as in the case of T^^, 

(Vj/Tq,^) = d^Iva^ - T^^Ixal3 - Tq^/i/A^ " ^^^luaX ■ (2-13) 

and in general requires the calculation of 1920 terms, of which fortunately many vanish by symmetry 
due to specific choice of metric. 

To start the process, one assumes for example that T^u has a perfect fluid form, 

T^,u = [p{t) + p{t) ]u^,u^ + g^^p{t) (2.14) 

for which one needs to compute the action of on T^j^, and then analytically continues the answer 
to negative fractional values of n = — l/2z/. Even in the simplest case, with (?(□) acting on a scalar 
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such as the trace of the energy-momentum tensor T^'^, one finds for the choice p{t) = pot^ and 
a{t) = ao the rather unwieldy expression 



n+l 



r(4 



Pot 



/3-2n 



(2.15) 



r(f + i-n)r(i^-n) 

with an integer n later analytically continued to n ^ — with v = ^. 

A more general calculation shows that a non-vanishing pressure contribution is generated in 
the effective field equations, even if one initially assumes a pressureless fluid, p{t) = 0. After a 
somewhat lengthy derivation one obtains for a universe filled with non-relativistic matter (p=0) 
the following set of effective Priedmann equations, 



k 



+ 



a' it) 
o?{t) ' o?{t) 

for the tt field equation, and 

k a?{t) 2d{t) 



Stt G(t) , , A 



l + ct(t/0'/" + 



pit) + 



Stt Go 



Q (t/*0 



+ ... 



p(t) + A 



(2.16) 



(2.17) 



a\t) a\t) a{t) 3 

for the rr field equation. In the above expressions the running of G appropriate for the RW metric 
is 



G{t) = Go 1 + 



5G{t) \ 



Gn 



l + ct 



t \ 1/'^ 

To, 



(2.18) 



Go J 

with Q of the same order as cq in Eq. (2.5), and to = C [13]; in the quoted reference it was estimated 
Q = 0.450 Co for the tensor box operator. Note that it is the running of G that induces an effective 
pressure term in the second (rr) equation, corresponding to the presence of a relativistic fluid due to 
the vacuum polarization contribution. One important feature of the new equations is an additional 
power-law acceleration contribution, on top of the standard one due to A. 



2.2 Introduction of the Wvac Peiranieter 



It was noted in [13] that the field equations with a running G, Eqs. (2.16) and (2.17), can 
be recast in an equivalent, but slightly more appealing, form by defining a vacuum polarization 
pressure pvac and density pvac, such that for the FLRW background one has 

Pvo.{t) = ^ p{t) p^acit) = Y-^ Pit) . (2.19) 

(jro O CrO 
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Consequently the source term in the tt field equation can be regarded as a combination of two 
density terms 

(l + ^) Pit) ^ Pit) + Pvacit) , (2.20) 

while the rr equation involves the new vacuum polarization pressure term 

l^-^pit)^Pvacit). (2.21) 

Form this viewpoint, the inclusion of a vacuum polarization contributions in the FLRW framework 
seems to amount to a replacement 

p{t) ^ p{t) + p^acit) Pit) ^ Pit) + p^acit) (2.22) 

in the original field equations. Just as one introduces the parameter w, describing the matter 
equation of state, 

pit) = w pit) (2.23) 

with w = for non-relativistic for matter, one can do the same for the remaining contribution by 
setting 

it) ■ (2.24) 

Then in terms of the two w parameters 

(w + Wyac^^^ Pit) = Pit) + Pvacit) (2.25) 

with, according to Eqs. (2.16) and (2.17) and following the results of [13], w^ac = 5 in a FLRW 
background. We should remark here that the calculations of [13] also indicate that w^ac = 5 is 
obtained generally for the given class of G(n) considered, and is not tied therefore to a specific 
choice of u, such as = |. 

The previous, slightly more compact, notation allows one to re- write the field equations for 
the FLRW background in an equivalent form, which we will describe below. First we note though 
that in the following we will restrict our attention mainly to a spatially flat geometry, k = 0. 
Furthermore, when dealing with density perturbations we will have to distinguish between the 
background, which will involve a background pressure ip) and background density ip), from the 
corresponding perturbations which will be denoted here by Sp and 6p. Then with this notation and 
for constant Go, the FLRW field equations for the background are written as 



a^{t) 

^ + 2^ = -8nGopit) + X. (2.26) 
a'^yt) a{t) 
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Now in the presence of a running (-?(□), and in accordance with the results of Eqs. (2.16) and 
(2.17), the modified FLRW equations for the background read 



-T77^ + ^^ = -SttGo iw + w^a^—^) p{t) + X, (2.27) 
a^{t) a[t) \ (jo / 

using the definitions in Eqs. (2.23) and (2.24), here with pvac{t) = WvacPvac{t)- 

We note here that the procedure of defining a Pvac and a Pvac contribution, arising entirely from 
quantum vacuum polarization effects, is not necessarily restricted to the FLRW background metric 
case [13]. In general one can decompose the full source term in the effective nonlocal field equations 
of Eq. (2.3), making use of 

as two contributions, 

i- G{U) T,, = (l + T,, = T,, + . (2.29) 

The latter involves the nonlocal part ^ 

T^T^^^T;.- (2.30) 
In addition, consistency of the full nonlocal field equations requires that the sum be conserved, 

(t^, + r;r) = . (2.31) 

It is important to note at this stage that the nature of the covariant d'Alembertian □ = gi^^ V^uVi^ 
is such that the result depends on the type of the object it acts on. Here T^i, is a second rank 
tensor (as in Eq. (2.13)), which causes a re-shuffling of components in T^i/ due to the matrix nature 
of both tensor □ and tensor G(n), and eventually accounts for the generation of a non- vanishing 
induced pressure term. This is clearly seen in the effective field equations of Eqs. (2.16) and (2.17), 
and in the ensuing definitions of Eq. (2.19). 

In general though one cannot expect that the contribution T^^"^ will always be expressible in 
the perfect fluid form of Eq. (2.14), even if the original T^u for matter (or radiation) has such a 
form. The former will in general contain, for example, non-vanishing shear stress contributions, 
even if they were originally absent in the matter part. Nevertheless the interesting question arises 
of whether, for example, Wyac = 5 continues to hold beyond the FLRW case treated above. In part 
this question will be answered affirmatively below, in the case of matter density perturbations. 

^Onc normally docs not include the l.h.s. field equation contribution +\g^ry as part of the r.h.s. matter part T^"'^, 
although it might be sensible to do so, given its large radiative (quantum) content [12]. We note here that the former 
is expected to contain the fundamental length scale ^ as well, in the form ~ -l-(3/^^) gi^i,. 



11 



3 Relativistic Treatment of Matter Density Perturbations 



Besides the modified cosmic scale factor evolution just discussed, the running of G(n) given in 
Eq. (2.28) also affects the nature of matter density perturbations on very large scales. In computing 
these effects, it is customary to introduce a perturbed metric of the form 

dr^ = dt^ - {Sij + hij) dx'dx^ , (3.1) 

with a{t) the unperturbed scale factor and hij{x, t) a small metric perturbation, and /too = hio = 
by choice of coordinates. As will become clear later, we will mostly be concerned here with the trace 
mode ha = h, which determines the nature of matter density perturbations. After decomposing the 
matter fields into background and fluctuation contribution, p = p + 6p, p = p + Sp, and v = v + (5v, 
it is customary in these treatments to expand the density, pressure and metric trace perturbation 
modes in spatial Fourier modes, 

Mx,i) = 5/9q(0e*'*-" 5p(x,0 = 5pq(i)e*''-" 

(5v(x,t) = (5vq(t)e^'i-" hij{:ic,t) = h^ij{t)e'''-'' (3.2) 

with q the comoving wavcmimbcr. Once the Fourier coefficients have been determined, the original 
perturbations can later be obtained from 

Spi^,t) = J ^0j^e-'^--Sp^{t) (3.3) 

and similarly for the other fluctuation components. Then the field equations with a constant Go 
(Eq. (2.2)) are given to zeroth order in the perturbations by Eq. (2.26), which fixes the three 
background fields a{t), p{t) and p{t). Note that in a comoving frame the four- velocity appearing 
in Eq. (2.14) has components = 1, = 0. To first order in the perturbations and in the limit 
q — > the field equations give 

a(t) 



ait) 



h{t) = 8TrGop{t)6{t) 



hit) + 3^ hit) = -2A7rGowpit)6it) (3.4) 
ait) 

with the matter density contrast defined as Sit) = 6 pit) /pit), hit) = huit) the trace part of hij, 
and It; = for non-relativistic matter. When combined together, these last two equations then 
yield a single equation for the trace of the metric perturbation, 

/i(i) + 2 ^ /i(i) = -87rGo(l + 3w)p(i)(^(t) • (3.5) 
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From first order energy conservation one has — ^ (1 + w) h{t) = d{t), whicli tlien allows one to 
eliminate h{t) in favor of d{t). This finally gives a single second order equation for the density 
contrast d{t), 

6{t) + 2 ^ 6{t) -AttG p{t) 5{t) = . (3.6) 

In the case of a running G(n) these equations need to be re-derived from the effective covariant 
field equations of Eq. (2.3), and lead to several additional terms not present at the classical level. 

Not surprisingly, as we shall see below, the correct field equations with a running G are not given 
simply by a naive replacement G G{t), which would lead to incorrect results, and violate general 
covariance. 



3.1 Zeroth Order Energy-Momentum Conservation 



As a first step in computing the effects of density matter perturbations one needs to examine 
the consequences of energy and momentum conservation, to zeroth and first order in the relevant 
perturbations. If one takes the covariant divergence of the field equations in Eq. (2.3), the left 
hand side has to vanish identically because of the Bianchi identity. The right hand side then gives 
+ T^y'^^ = 0, where the fields in T^"'^ can be expressed, at least to lowest order, in terms 
of the pyac and p^ac fields defined in Eqs. (2.19) and (2.24). The first equation one obtains is the 
zeroth (in the fluctuations) order energy conservation in the presence of (?(□), which reads 



a{t) 



ait) 

For w = and Wvac 



{l + w) + + 



5G{t) 



,M,^,M,(l,^|p)^(,).0. (3.7, 



I this reduces to 



, a(t) ^ ^ a{t) Sa(t) ^ SG(t) 



m = , 



ait) a{t) Go ' Go 
or equivalently in terms of the variable a{t) only 

p-(a) + (l + ^)p»=0. 



(3.8) 



3 4 6G{a) 6G'ia) 
a a Go Go 



(3.9) 



In the absence of a running G these equations reduce to the ordinary mass conservation equation 
for w = 0, 

m = -^^m- (3.10) 



ait) 
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It will be convenient in the following to solve the energy conservation equation not for p{t), but 
instead for p{a). This requires that, instead of using the expression for G{t) in Eq. (2.18), one uses 
the equivalent expression for G{a) 

C,„) = 0.(l.^|5)). „,,^fl...(;^)%.... (3.11) 

In this last expression the power is = 3/2i/, since from Eq. (2.18) one has for non-relativistic 
matter a(t)/aQ ^ {t/tof/^ in the absence of a running G. In the following we will almost exclusively 
consider the case = | [8] for which therefore 7,^ = 9/2. ^ Then in the above expression Cq q if 
ag is identified with a scale factor appropriate for a universe of size ^; to a good approximation this 
should correspond to the universe "today" , with the relative scale factor customarily normalized at 
such a time to a/a^ = 1. Consequently, and with the above proviso, the constant Ca in Eq. (3.11) 
can safely be taken to be of the same order as the constant cq appearing in the original expressions 
for G{U) in Eq. (2.28). 

Then the solution to Eq. (3.8) can be written as 

{ f da SG{a) SG'{a)\] 

p{a) = const, exp |- y — (^3 + + a | , (3-12) 

or, more explicitly, as 

Vl + ^«te) ; l + (l + 7.-^)c„(^) 

with p{a) normalized so that p{a = ao) = po- For Cq = the above expression reduces of course to 
the usual result for non-relativistic matter, 

m=Po (^)' . (3.14) 

Furthermore, here one also finds that the zeroth order momentum conservation equation is identi- 
cally satisfied, just as in the case of constant G. 



3.2 Zeroth Order Field Equations with Running G(n) 



''This implicitly assumes that the cosmic evolution is largely matter dominated, if p = wp then a{t)/ao = 
{t/tof'^^''+'"\ In the opposite regime where a cosmological constant can eventually prevail one has instead 
a{t)/ao = exp^/\/^{t-to). Then j- = I + j-^/\\og f - and for to ^ C and ^J\ ~ (, one has simply ^ = l-|-log ^. 
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The zeroth order field equations with the running of G included were already given in Eq. (2.27). 
One can subtract the two equations from each other to get an equation that does not contain A, 



a\t) a{t) 



{l + w) + {l + 



SGit) 



-p(t) . (3.15) 



Go 

Alternatively, from Eqs. (2.27) one can obtain a single equation that only involves the acceleration 
term with a{t), 

n(i\ r ^n(f^^ 

p{t) + A . (3.16) 



3^ = -4.Go 

a{t) 



{l + 3w) + {l + W^ac) 



Go 

It is also rather easy to check the overall consistency of the energy conservation equation, Eq. (3.8), 
and of the two field equations in Eq. (2.27). This is done by (i) taking the time derivative of the 
first tt equation in Eq. (2.27), (ii) replacing terms involving p hy p using the energy conservation 
equation, Eq. (3.8), and (iii) finally by substituting again the result of the first (tt) equation into 
Eq. (2.27) to obtain the second (rr) equation in Eq. (2.27). 



3.3 Effective Energy-Momentum Tensor p^ac? Pvac 

The next step consists in obtaining the equations which govern the effects of small field per- 
turbations. These equations will involve, apart from the metric perturbation hij, the matter and 
vacuum polarization contributions. The latter arise from (see Eq. (2.29)) 

(l + T,, = T,, + (3.17) 

with a nonlocal T^"'^ = {5G{U)/Gq) T^^. Fortunately to zeroth order in the fluctuations the results 
of Ref. [13] indicated that the modifications from the nonlocal vacuum polarization term could 
simply be accounted for by the substitution 

p{t) p{t) + p^ac{t) Pit) ^ Pit) + Pyacit) . (3.18) 

Here we will apply this last result to the small field fluctuations as well, and set 

Spqit) ^ 6pq{t) +6pqyac{t) 6pq{t) ^ 6pq{t) + dp^^acit) ■ (3.19) 

The underlying assumptions is of course that the equation of state for the vacuum fluid still remains 

roughly correct when a small perturbation is added. Furthermore, just like we had p{t) = w p{t) 
(Eq. (2.23)) aad pyac{t) = Wyac Pvac{t) (Eq. (2.24)) with Wyac = |) we now write for the fluctuations 

Spdit) =w6pci{t) Spqvacit) = Wyac Spqvacit) , (3.20) 
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at least to leading order in the long wavelength limit, q ^ 0. In this limit we then have simply 

5p{t) = W dp{t) SpvaS) = W^ac ^PvaS) = Wvac —^Sp{t) , (3.21) 

with G{t) given in Eq. (2.18), and we have used Eq. (2.19), now appUed to the fluctuation 6pvac{t), 

Spvacit) = ^-^Sp{t) + ... (3.22) 

where the dots indicate possible additional 0{h) contributions. 

A bit of thought reveals that the above treatment is incomplete, since G{0) in the effective field 

equation of Eq. (2.3) contains, for the perturbed RW metric of Eq. (3.1), terms of order hij, which 
need to be accounted for in the effective T^^^c- Consequently the covariant d'Alembertian has to be 
Taylor expanded in the small field perturbation hij, 

a{g) = + □(i)(/i) + 0{h'^) , (3.23) 

and similarly for G{n) 



G{a) = Go 



CO / 1 



1 + 0(0) + □(!) (^) + 0(^2) ) + ■ ■ ■ 



(3.24) 



which requires the use of the binomial expansion for the operator (A+B)^^ = A^^ — A^^ B + . . .. 
Thus for sufficiently small perturbations it should be adequate to expand (?(□) entering the effective 
field equations in powers of the metric perturbation hij . Since a number of subtleties arise in this 
expansion, we shall first consider the simpler case of a scalar box, where some of the issues we 
think can be clearly identified, and addressed. After that, we will consider the more complicated 
case of the tensor box. This will be followed by a determination of the effects of the running of G 
on the relevant matter and metric perturbations, by the use of the modified field equations, now 
expanded to first order in the perturbations. 



3.4 0(h) Correction using Scalcir Box 



In this section the term 0{h) in Sp-^ac of Eq. (3.19) will be determined, using a set of formal 
manipulations involving the covariant scalar box operator. Instead of considering the full field 
equations with a running G(n), as given in Eq. (2.3), 



^ + A = Stt Go (l + T^u (3.25) 
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we will consider here instead the action of a scalar on the trace of the field equations for 

A = 0, 



R 



-SttGo 1 + 



5G{p, 
Go 



or equivalently, by having the operator (?(□) act on the left hand side, 



Gn 



For a perfect fluid one has simply T^' 



-p, which then gives [13] 



Go 



or equivalently 



Gn 



i.cll) .... 



p{t) = G{t)p{t) 



(3.26) 



(3.27) 



(3.28) 



(3.29) 



with ct — 0.785 Co, and to = ? [13] (in the tensor box case a slightly smaller value was found, 
Ct 0.450 Co). The two terms in Eq. (3.29) are of course recognized, up to a factor of Go, as the 
combination 

p{t) + Pvacit) (3.30) 

of Eq. (3.18), with pyac{t) = SG{t)/Go ■ p{t). Thus the zeroth order result obtained by the use of 
the scalar d'Alembertian acting on the trace of the field equations is consistent with what has been 

used so far for G[t). 

To compute the higher order terms in the hifs appearing in the metric of Eq. (3.1) one needs 
to expand G(n) according to Eq. (3.24) giving 



G(n) = Go 



Here we are interested in the correction of order /ijj, when the above operator acts on the scalar 
-p. This would then give the correction 0{h) to Sp^ac^ namely the second term in 



A 

-'a 



, , 5G(a(°)) , , , 5G(u)(h) , , 



(3.32) 



Go Go 

with the first term being simply given in the FLRW background by 6G{t)/Go -5 pit). Here the 0{h) 
correction is given explicitly by the expression 



5G{n){h) 1 CO 1 / 1 \V2'^ 



Go 



P = -TT- 



(3.33) 



The effect of the (U^^^) ^1"^^ term is essentially to make the coupling time dependent, i.e. to 
correctly reproduce the required overall time-dependent factor 5G(t)/GQ. 
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Now the scalar d'Alembertian □ = g^'^ 'V^'Vu acting on scalar functions S{x) has the form 

□ S{x) = ^df, g^'^d. Six) (3.34) 

In the absence of hij fluctuations this gives for the metric in Eq. (3.1) 

□ (o)5(x) = 4tV25-3-5-5 ^ (-df -3-dt) S(t) , (3.35) 
a^^ a \ a J 

where in the second expression we have used the properties of the RW background metric: we only 

need to consider functions that are time dependent, so that S{x,t) S{t). To first order in the 

field fluctuation hij of Eq. (3.1) one computes 



□«(/t) Six) = --hs- -K^ dls + ^ i-d^ K^) ■d^s + -^ d^h ■d^s + ... (3.36) 



with the trace hit) = hxxit) + hyyit) + hzzii)- But for a function of the time only one obtains 

n(^)ih)pit) = ~hit)Sit). (3.37) 

Thus to first order in the fluctuations one obtains the expression 

^ • □W(^) . iSGp) = Z^^^ • ^ (3 ^ - ^g) -p (3.38) 

where use has been made of the zeroth order mass conservation equation in Eq. (3.10). Note that 
this result also correctly incorporates the effect of G(n('^)) on functions of t, as given for example 
in Eq. (3.28), which ensures the proper running of (5G(t). 

Now in our treatment we are generally interested in mass density and metric perturbations 
around a near-static background described by a/ a = Hit), and pit). For these we expect the 
relevant time variations in Sp and h to be somewhat larger than for the background itself. Thus 
for sufficiently slowly varying background fields we retain only /i(t) and its derivatives, and for a 
sufficiently slowly varying /i(t) only hit) and the lowest derivatives. Then the factors of a/ a are 
seen to cancel out at leading order between numerator and denominator in Eq. (3.38), and one is 
left simply with 

^ □(^^ ih) ■ 6Git) Pit) = - I 6Git) hit) Pit) + ... (3.39) 



□ (0) ^ ^ \ jy\ J 2 
Putting everything together, one finds for the O(^) correction 

SGia)ih) , , 1 bGit) , , , , , , , , 

The scalar box calculation just described allows one to compute the correction 0(/t) to Spvacit) in 
Eq. (3.32), and leads to the following 0(/i) modification of Eq. (3.22) 

Spvacit) = ^ 6pit) + ±cj-^ hit) Pit) (3.41) 
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and similarly from 5pvac{t) = w^ac^ Pvadt) , 

Spvacit) = W,ac Spit) + ^ ^ ^t) (3.42) 

with Wyac = I- The second 0(/i) terms in both expressions account for the feedback of the metric 
fluctuations h on the vaccum density 5p„ac and pressure Sp^ac fluctuations. 

The potential flaw with the preceding argument is that it assumes that certain very specific 
functions of the background stay constant, or at least very slowly varying. In the case at hand 
this was a/a = i?(a) const, and p ^ const., which in principle is not the only possibility, and 
would seem therefore a bit restrictive. A slightly more general approach, and a check, to the 
evaluation of the expression in Eq. (3.38) goes as follows. One assumes instead a harmonic time 
dependence for the metric fluctuation hit) = hg e*'^*, and similarly for a(t) = e*^*, pit) = pQ e*^*, 
and 6Git) = SGo e*^*; different frequencies for a and p could be considered as well, but here we will 
just stick with the simplest possibility. Then from the last expression in Eq. (3.38) one has 

lh(3^5G-S-G)p= , , ^^^, -i-rudGhp) . (3.43) 



In the limit w 3> F, corresponding to h/h ^ a/ a, one obtains for the above expression 

- - SGit) hit) pit) - f-i) SGit) hit) pit) , (3.44) 



00 \a hy 

after substituting back h/h = icj and a/a = iV in the last expression. Then Spvacit) in Eq. (3.41) 
now involves the quantity Ch 

Ch = -i. (3.45) 

a h 

At first this last factor (a function and not a constant) would seem rather hard to evaluate, and 
perhaps not even close to constant in time. But a bit of thought reveals that, to the order we are 
working, one can write 

ha^ dloghja) ^ d log 6ia) ^ 

h a 5 log a 9 log a ' 

where (5(a) is the matter density contrast, and /(a) the known density growth index [17]. In the 

absence of a running G (which is all that is needed, to the order one is working here) an explicit 

form for /(a) is known in terms of derivatives of a Gauss hypergeometric function, which will 

be given below. One can then either include the explicit form for /(a) in the above formula for 

Spvacit)-, or use the fact that for a scale factor referring to "today" a/ao ^ 1, and for a matter 

fraction CI ^ 0.25, one knows that /(a = oq) — 0.4625, and thus in Eq. (3.41) one obtains the 

improved result Ch — 2.1621. This can then be compared to the earlier result, which gave Ch — 1/2. 
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A similar analysis can now be done in the opposite, but in our opinion less physical, w <C T 
limit, for which one now obtains for the expression in Eq. (3.43) 

This new limit is less physical because of the fact that now the background is assumed to be varying 
more rapidly in time than the metric perturbation itself, d/a S> h/h. For dpyac{i) one then obtains 
a similar expression to the one in Eq. (3.41), with a different coefficient 

1 a h ,„ 
c/. = TTT V -r 3.48 
10 a h 

still involving the quantity {a/a){h/h) = f(a). By the same chain of arguments used in the previous 
paragraph one can now either include the explicit form for /(a) in the formula for 5pvac{t)^ or: use the 
fact that for a scale factor referring to "today" a/oo ~ 1 and a matter fraction $7 0.25 one knows 
that /(a = ao) 2± 0.4625. In this case one then has in Eq. (3.41) ~ (1/10) x 0.4625 = 0.0463. 
One disturbing, but not entirely surprising, general aspect of the whole calculation in this second 
a; ^ r limit (as opposed to the previous treatment in the opposite limit) is its rather significant 
sensitivity, in the final result, to the set of assumptions initially made about the time development 
of the background as specified by the functions a{t) and p{t). Therefore in the following we shall 
not consider this limit further. 

To summarize, the results for a scalar box and a slowly varying background, h/h S> d/a, give the 
0{h) corrected expression for 6pvac{t) in Eq. (3.41) and 6pvac{t) = Wyac Spvac{t), with c/i ~ +2.1621. 



3.5 0(h) Correction using Tensor Box 



The results of Eqs. (3.32) and (3.41) for the vacuum polarization contribution, 

6pvac{t) = ^ Sp{t) + J- c, ^ h{t) -p{t) (3.49) 

and similarly for 5pvac{t) = Wvac ^Pvac{t) with w^ac = \^ were obtained using a scalar d'Alembertian 
to implement G(n) by considering the trace of the field equation, Eq. (3.26). In this section wc will 
discuss instead the result for the full tensor d'Alembertian, as it appears originally in the effective 
field equations of Eqs. (2.3) and (3.25). 

Now the d'Alembertian operator □ = g^^V acts on the second rank tensor T^j, as in 
Eq. (2.13), and should therefore be regarded as a four by four matrix, transforming T^^, into 
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[□T]^,^. Indeed it is precisely this matrix nature of □, and therefore of that accounts for the 

fact that a vacuum pressure is induced in the first place, leading to a Wvac 7^ 0. 

To compute the correction of 0(h) to 5pvac{t) one needs to consider the relevant term in the 
expansion of (1 + 8G(p) j Gq) Tf^v, which we write as 



-^J-. ^^^-T - -J- □(!)(/,). r^V^''-r (3 50) 

2uum Go ~ 21/ ^V- 0(0) ^ [nWj ^^'-^ ^'^■^^1 

This last form allows us to use the results obtained previously for the FLRW case in [13], namely 

-^^yr^T,, = (3.51) 

with here 

T^^T = \Pvac{t)+Pvac{t)\u^U^ + g^,^p^ac{t) (3.52) 

and (see Eqs. (2.19) and (2.30)), to zeroth order in h, 

, . 6G(t) _, . , , SG(t) _, . , , 

PvoAt) = pit) Pvac[t) = Wyac P{t) ■ (3.53) 

with Wvac = |. Therefore, in light of the results of Rcf. [13], the problem has been dramatically 
reduced to just computing the much more tractable expression 

-2^BM-°^'^(^)-^mT' (3.54) 
and in fact the only ordering for which the expression ((^G(n)/Go) T^j^ is calculable within rea- 
sonable effort. Still, in general the resulting expression for • □(^)(/i) is rather complicated if 
evaluated for arbitrary functions, although it does have a structure similar to the one found for the 
scalar box in Eq. (3.38). 

Here we will resort, for lack of better insights, to a treatment that parallels what was done 
before for the scalar box, where one assumed a harmonic time dependence for the metric trace 
fluctuation h{t) = h^e^*', and similarly for a{t) = ao e*'"* and p{t) = poe*'"*. In the limit a; S> F, 
corresponding to h/h 2> a/a, one finds for the fiuctuation Sp^acit) in Eq. (3.41) 

5pvac{t) = ^ Spit) + ^Ch ^ hit) Pit) . (3.55) 

The 0(/i) correction factor c/j for the tensor box is now given by 

Ch = — - J , (3.56) 

o a h 

with all other off-diagonal matrix elements vanishing. Furthermore one finds to this order, but only 
for the specific choice w^ac = 5 the zeroth order TJ^^'^, 

Spvacit) = ^Spvacit) (3.57) 
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i.e. the 0{h) correction preserves the original resuh w^ac = 5- other words, the first order result 
0(h) just obtained for the tensor box would have been somewhat inconsistent with the zeroth order 
result, unless one had w^ac = g to start with. Now, one would not necessarily expect that the first 
order correction could be still be cast in the form of the same equation of state pvac — ^Pvac as the 
the zeroth order result, but it would nevertheless seem attractive that such a simple relationship 
can be preserved beyond the lowest order. 

As far as the magnitude of the correction Ch in Eq. (3.56) one can argue again, as was done in 
the scalar box case, that from Eq. (3.46) one can relate the combination {h/h){a/a) to the growth 
index /(a). Then, in the absence of a running G (which is all that is needed here, to the order 
one is working), an explicit form for /(a) is known in terms of suitable derivatives of a Gauss 
hypergeometric function. These can then be inserted into Eq. (3.56). Alternatively, one can make 
use again of the fact that for a scale factor referring to "today" a/ao ~ 1, and for a matter fraction 
n « 0.25, one knows that /(a = uq) ^ 0.4625, and thus in Eq. (3.41) Ch (11/3) x 2.1621 = +7.927. 
This last result can then be compared to the earlier scalar result which gave c/j 2± +2.162 using 
the same set of approximations (slowly varying background fields) . It is encouraging that the new 
correction is a bit larger but not too different from what was found before in the scalar box case. 
Note that so far the sign of the 0{h) correction is the same in all physically relevant cases examined. 

Next, as in the scalar box case, one can do the same analysis in the opposite, but less physical, 
limit w <C r or <C a/a. One now obtains from the tt matrix element the 0{h) correction in 
the expression for dp^ac given in Eq. (3.41), namely 

^cj-^ hit) Pit) ■ (3.58) 



with a coefficient 



121 a;2 121 /a\2 ^ 



Similarly for the ii matrix element of the 0(/i) correction one finds 



with 



^a\t)d^^-^hit)-pit). (3.60) 



4 = (3.61) 



giving now the Spmcih) correction. Again all off-diagonal matrix elements are equal to zero. It 
seems therefore that this limit, lo <^T or h/h <^ a/ a, leads to rather different results compared to 
what had been obtained before: the only surviving contribution to 0(/i) is a rather large pressure 
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contribution, with a sign that is opposite to all other cases encountered previously. Furthermore 
here the relationship Wvac = | is no longer preserved to 0{h). But, as emphasized in the previous 
discussion of the scalar box case, this second limit is in our opinion less physical, because of the 
fact that now the background is assumed to be varying more rapidly in time than the metric 
perturbation itself, d/a S> h/h. Furthermore, as in the scalar box calculation, one disturbing 
but not entirely surprising general aspect of the whole calculation in this second a; ^ F limit, 
is its extreme sensitivity as far as magnitudes and signs of the results are concerned, to the set 
of assumptions initially made about the time development of the background. As a final sample 
calculation let us mention here the case, similar to what was done originally for the scalar box, where 
one assumes instead a/a = H{a) « const, and p « const., which, as we mentioned previously, seems 
now a bit restrictive. Nevertheless we find it instructive to show how sensitive the calculations are 
to the nature of the background, and in particular its assumed time dependence. In the notation 
of Eqs. (3.58), (3.59) and (3.61) one finds in this case 

_ 625 u;2 _ 625 1 /i , _ 4 

Again here the pressure contribution Sp^acih) is the dominant contribution, the 5pyac{h) part being 
negligible, 0{h). For the reasons mentioned, in the following we will no longer consider this limit 
of rapid background fluctuations any further. 

To summarize, the results for a scalar box and for a very slowly varying background, h/h^ a/ a, 
give the 0{h) corrected expression for 6p-uac{t) in Eq- (3-41) and Spyac{t) = w^ac^Pvacit) with 
Cfi 2± +2.162, while the tensor box calculation, under essentially the same assumptions, gives the 
somewhat larger result ~ +7.927. From now on, these will be the only two choices we shall 
consider here. 

3.6 First Order Energy-Momentum Conservation 

The next step in the analysis involves the derivation of the energy-momentum conservation to 
first order in the fluctuations, and a derivation of the relevant field equations to the same order. 
After that, energy conservation will be used to eliminate the h field entirely, and thus obtain a 
single equation for the matter density fluctuation 5. 

The results so far can be summarized as follows. For the metric in Eq. (3.1), and in the limit 
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q 0, the field equations in Eq. (2.3) can now be written as as 

R^^. - 5 ffM- ^ + AffM- = Stt Go (t^. + r;^) , (3.63) 

with T^j"*^ = {5G(0)/Gq)Thi,. Here T^^ describes the ordinary matter contribution, in the form 
of a perfect fluid as given in Eq. (2.14), here with p = wp and it; 2± 0, while T^^*^ describes the 
additional vacuum polarization contribution 

^/^r = [Pvac{t) + Pvac{t) ]u^U^ + g^,uPvac{t) (3.64) 

with pvac = WvacPvac and Wvac = 5, as in Eq. (2.24). Furthermore, each field now contains both a 
background and a perturbation contribution, 

p{t) = p{t) + 5p{t) p{t) = w p{t) (3.65) 

and similarly 

Pvac\J') Pvac (t) + 5pyac{t) 

PvacV^j '^vac Pvac 

(t) . (3.66) 

Prom Eq. (2.19) one has 



Pvac{t)=^-^p{t), (3.67) 



while from Eq. (3.41) on has 



6pvac{t) = ^ Sp{t) + J- c, ^ h{t) Pit) (3.68) 

and similarly Spvacif) = Wvac ^ Pvac{t) . The second 0{h) terms in both expressions physically 
account for the feedback of the metric fluctuations h on the vaccum density Sp^ac and pressure 
Spvac fluctuations. In light of the discussion of the previous section, wc will limit our derivations 
below to the case of constant c/j; the case of a non-constant as in Eq. (3.46) can be dealt with 
as well, but the resulting equations are found to be quite a bit more complicated to write down. 

Consequently all quantities in the effective field equations of Eq. (3.63) have been specified to 
the required order in the field perturbation expansion. First we will look here at the implications 
of energy-momentum conservation, (t^i, + T^"*^^ = 0, to first order in the fluctuations. The 
zeroth order energy conservation equation was already obtained in Eq. (3.7), and its explicit solution 
for p{a) given in Eq. (3.13). After defining the matter density contrast S{t) as the ratio S{t) = 
Sp{t)/p{t), the energy conservation equation to first order in the perturbations is found to be 



1 A, , , (5G(t)\ 1 5G{t) ■ 

--((l + .) + (l + ....) ^j-^c.^ 

+ ^^Ch y6[W W,ac) ^^^^ j 
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h{t) 

6G{t) 



1 + 



Go 



8{t) . (3.69) 



In the absence of a running G {6G{t) = 0) this reduces simply to — ^ (1 + w) h{t) = S(t), and thus 
to the standard result for the metric trace perturbation in terms of the density contrast 



(1 + ^) h{t) = S{t) 



(3.70) 



This last result then allows us to solve explicitly, at the given order, i.e. to first order in the 
fluctuations and to first order in SO, for the metric perturbation h{t) in terms of the matter 
density fluctuation S{t) and 5{t), 

2 



h{t) 



1 + w 
1 



1 A , 1 \ SG(t) 

1 + w V 2z/, 



2u {l + wY 



3 {W- Wyac) 



Go 

a{t) SG{t) SG{t) 



a{t) Go 



Go 



m 



(3.71) 



Similarly, by differentiating the above relationship, an expression for h{t) in terms of S and its 
derivatives can be obtained as well. 



3.7 First Order Field Equations 



To first order in the perturbations, the tt and ii effective field equations become, respectively, 
^ h{t) - Stt Go ^ c, ^ m h{t) = 8nGo(l + ^) m 6{t) (3.72) 

and 

Mt) + 3 4S hit) + 247r Go ^ w,ac ^ m h{t) = - 247r Go (w + w,ac ^) p(t) S{t) 
a[t) zu Ctq \ (^0 J 

(3.73) 

In the second ii equation, the zeroth order ii field equation of Eq. (2.27) has been used to achieve 
some simplification. 

As a final exercise, it is easy to check the overall consistency of the first order energy conservation 
equation of Eq. (3.69), and of the two field equations given in Eqs. (3.72) and (3.73). To do so, one 
needs to (i) take the time derivative of the tt equation in Eq. (3.72); (ii) get rid of p consistently 
by using energy conservation to zeroth order in 6G and in the fluctuations from Eq. (3.69) for 
terms of order 6G times a fluctuation, combined with the use of energy conservation to first order 
in SG but without fluctuations as in Eq. (3.8) for the terms that are already of first order in the 
fluctuations; (iii) eliminate the 5 terms using the energy conservation equation to first order in SG 
without field fluctuations (Eq. (3.8)) for terms proportional to SG times a fluctuation, and using 
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the energy conservation equation to first order in 6G and in the fluctuation (again Eq. (3.69)) for 
terms of zeroth order in the fluctuations; (iv) use the combination of Eqs. (2.27) that does not 
contain A, Eq. (3.15), to get rid of aja terms; (v) Finally use the tt equation for the fluctuation, 
Eq. (3.72), to eUminate some terms proportional to p times a fluctuation so as to finally obtain the 
second ii field equation Eq. (3.73). 



3.8 Matter Density Contrast Equation in t 



To obtain an equation for the matter density contrast 5{t) = 5p{t)/p{t) one needs to eliminate 
the metric trace field h{t) from the field equations. This is first done by taking a suitable linear 
combination of the two field equations in Eqs. (3.72) and (3.73), to get the equivalent equation 
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Go 



■p{t)h{t) 



■SttGo 



(l + 3w;) + (l + 3 



5G{t) 



Go 



mm- 



(3.74) 



Then the first order energy conservation equations to zeroth (Eq. (3.70)) and first (Eq. (3.71)) 
order in SG allow one to completely eliminate the h, h and h field in terms of the matter density 
perturbation 6{t) and its derivatives. The resulting equation reads, for -u; = and Wyac = ^, 



m 



+ 



- 47r Go ( 1 + 



2u 



1 5G{t) \ 
3 Go J" 
7 5G(t) 

3~Gr 

a^{t) 5G{t) 



2u 



a{t) 8G{t) ^ jG{t)\ 
Go ) 



Y.-^""'- Go 



a{t) Go 
5G{t) 



m 



m 



^^a(t) 5G(t) _^a(t) SG{t) 



a^{t) Go a{t) Go a{t) Go 



SG{t) \ 
Go 



5{t) = 



(3.75) 



This last equation then describes matter density perturbations to linear order, taking into account 
the running of G(n), and is therefore the main result of this paper. The terms proportional to 
Cfi, which can be clearly identified in the above equation, describe the feedback of the metric 
fluctuations h on the vaccum density Sp^ac a^nd pressure Sp^ac fluctuations. The equation given 
above can now be compared with the corresponding, much simpler, equation obtained for constant 
G, i.e., for G — )■ Go and still w = (see for example [18] and [17]) 



S{t) + 2- 5{t) - 47r Go p{t) 6{t) = 



(3.76) 
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from which one obtains for the growing mode 

2/3 



S^{t) = S^{to) (^- j , (3.77) 



which is the standard result in the matter-dominated era. 



3.9 Matter Density Contrast Equation in a(t) 



It is common practice at this point to write an equation for the density contrast d{a) as a 
function not of t, but of the scale factor a{t). This is done by utilizing the following simple 
derivative identities 



f{t) = aH{a) 



f{t) = a^H\a) 



dlnH{a) 



+ 



1 



dfja) 

da 
df{a) 



d^f{a) 



(3.78) 
(3.79) 



da ' aj da ' ^ ^ da? 

where / is any function of t, and H = a{t)/a{t) the Hubble constant. This last quantity can be 
obtained from the zeroth order tt field equation 
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Often this last equation is written in terms of current density fractions, 

2 / i \ 2 
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n (i + zy + i^R {i + zY + n^ 



(3.80) 



(3.81) 



with a/ao = 1/(1 + 2;) where z is the red shift, and ao the scale factor "today". Then Ho is the 
Hubble constant evaluated today, Q, the (baryonic and dark) matter density, the space curvature 
contribution corresponding to a curvature k term, and Q\ the dark energy or cosmological constant 
part, all again measured today. In the absence of spatial curvature k = one has today 
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(3.82) 



In terms of the scale factor a{t) the equation for matter density perturbations for constant G = Gq, 
Eq. (3.76), becomes 
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The quantity H{a) is most simply obtained from the FLRW field equations 

H{a) = JfGop{a) + ^ , 



(3.83) 



(3.84) 
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with the matter density given in Eq. (3.14), which can in principle be solved for a(t), 

t-to= [ , ^" (3.85) 
J aVf Gop-o(^)' + | 

It is convenient at this stage to introduce a parameter describing the cosmological constant 
fraction as measured today, 

^=87g^ - - ~jr • ^^-^^^ 

While the following discussion will continue with some level of generality, in practice one is mostly 
interested in the observationally favored case of a current matter fraction O 0.25, for which 3. 
In terms of the parameter 9 the equation for the density contrast 5{a) for constant G can then be 
recast in the slightly simpler form 

dH{a) 3(l + 2a3g) d 6{a) _ 3 _ 

da^ ^2a{l + aH) da 2 (1 + ^) ^ ^ " ^ ' 

A general solution of the above equation is given by a linear combination of the two solutions 

Soia) = ci ■ + ^ a-^/^ + C2 ■ a ■ 2F1 Q,l;^;-a^e^ (3.88) 

where ci and C2 are arbitrary constants, and 2-^1 is the Gauss hypergeometric function. The 
subscript in So{a) is to remind us that this solution is appropriate for the case of constant 
G = Gq. Since one is only interested in the growing solution, the constant c\ = 0. 

To evaluate the correction to 5o(«) coming from the terms proportional to one sets 

5{a) oc (5o(a) [ 1 + c„ J^(a) ] , (3.89) 

and inserts the resulting expression in Eq. (3.75), written now as a differential equation in a(i), 
after using Eqs. (3.78) and (3.79) to replace 

a{t) = aH 

a,) . .^^^(i^.i). ,3.0) 

One only needs to determine the differential equations for density perturbations 6 up to first order 
in the fluctuations, so it will be sufficient to obtain an expression for Hubble constant H from the 
tt component of the effective field equation to zeroth order in the fiuctuations, namely the first of 
Eqs. (2.27). One has 

H(u) = Go(l + ^) m + \ (3.91) 
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with G{a) given in Eq. (3.11) and p{a) given in Eq. (3.12). ^ In this last expression the exponent 
is 7t/ = 3/2z^ ~ 9/2 for a matter dominated background universe, although more general choices, 
such as "f,y = 3(1 + w)l1v or even the use of Eq. (3.85), are possible and should be explored (see 
discussion later). Also, Ca ~ cj if ao is identified with a scale factor corresponding to a universe of 
size ^; to a good approximation this corresponds to the universe "today", with the relative scale 
factor customarily normalized at that time to a/oo = 1. In [13] it was found that in Eq. (2.18) 
ct 0.785 Co in the scalar box case, and ct 2± 0.450 cq in the tensor box case; in the following we 
will use the average of the two values. 

After the various substitutions and insertions have been performed, one obtains, after expanding 
to linear order in ao, a second order linear differential equation for the correction -F(a) to S{a), as 
defined in Eq. (3.89). Since this equation looks rather complicated for general 5G{a) it won't be 
recorded here, but it is easily obtained from Eq. (3.75) by a sequence of straightforward substitutions 
and expansions. The resulting equation can then be solved for J^{a), giving the desired density 
contrast (5(a) as a function of the parameter J7. 

Nevertheless with the specific choice for G{a) given in Eq. (3.11) an explicit form for the equation 
for (5(a) reads: 

^ + ^(<.)^+B(„)*(a) = 0. (3.92) 

with the two coefficients given by 
_ 3(l + 2a3g) _ 

c„(9a3(l+7^)^z/+a>(6c,,7^(l+27^)(l+a30)V(-9a3^+7^(l+a30)(3+27^(l+a3e)))i/)) 

6ai^7^. (l+a36')^ 

(3.93) 

and 

3 



c„(3a3(l+7^)gt/+aT'-(cft7^(2+7^)(l+a3g)(-l+27^+2a3(l+7^)g) + (47^+a3(-3+47^)g)i/)) 



2i/7^a2(l + a3e) 



(3.94) 



and the variable a considered just as a stand-in for what should really be the variable a/ao. To 

obtain an explicit solution to the above equation one needs to know the coefficient and the 



*We have noted before that Eq. (3.91) is suggestive of a small additional matter contribution, Q„ 
{8'K/3)SG{a)po/Hi, to the overall balance in Eq. (3.82). 



29 



exponent in Eq. (3.11), whose likely values are discussed above and right after the quoted 
expression for G{a). For the exponent u one has ~ ^, whereas for the value for Ch one finds, 
according to the discussion in the previous sections, c/j ~ 7.927 for the tensor box case. Furthermore 
one needs at some point to insert a value for the matter density fraction parameter 6 as given in 
Eq. (3.86), which based on current observation is close to = (1 — r2)/r2 2± 3. 



4 Relativist ic Growth Index with G(n) 



The solution of the above differential equation for the matter density contrast in the presence of 
a running Newton's constant G{D) leads to an explicit form for the function S{a) = So{a) [l+CaT{a)] . 
Prom it, an estimate of the size of the corrections coming from the new terms due to the running 
of G can be obtained. It is clear from the previous discussion, and the form of (?(□), that such 
corrections will become increasingly important in the present era t ~ to or a ~ cq. When discussing 
the growth of density perturbations in classical General Relativity it is customary at this point to 
introduce a scale-factor-dependent growth index f{a) defined as 

/(«)-^. (4.1) 

a ma 

which is in principle obtained from the differential equation for any scale factor a{t). Nevertheless, 
here one is mainly interested in the neighborhood of the present era, a{t) ^ qq. One therefore 
introduces today's growth index parameter 7 via 

a In (5(a) 



/(a = ao) = 



5 In a 

so that the exponent 7 itself is obtained via 

In/ 

7 = 



Inn 



= W , (4.2) 

a=aQ 



(4.3) 



a=ao 

The solution of the above differential equation for S{a) then determines an explicit value for the 
growth index 7 parameter, for any value of the current matter fraction ft. In the end, because of 
observational constraints, one is mostly interested in the range 0, ~ 0.25, so the following discussion 
will be limited to this case only, although from the original differential equation for 6 (a) one can 
in principle obtain a solution for any sensible $7. Numerically the differential equation for S{a) can 
in principle be solved for any value of the parameters. In practice we have found it convenient, 
and adequate, to obtain the solution as a power series in either f2 or 1 — In the first case the 
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resulting series is asymptotic and only slowly convergent around Q ~ 0.25, while in the latter case 
the convergence is much more rapid. In this last case we have carried therefore the expansion up to 
eighth order, which gives the answers given below (see also Figures 1-4) to an accuracy of several 
decimals. 

It is known that in the absence of a running Newton's constant G {G ^ Go, thus Ca = 0) one 
has /(a = ao) = 0.4625 and 7 = 0.5562 for the standard AC DM scenario with ^ = 0.25 [17]. On 
the other hand, when the running of G{n) is taken into account, one finds from the solution to 
Eq. (3.75) for the growth index parameter 7 at J7 = 0.25 the following set of results. 

For the tensor box case discussed in Sec. (3.5) one has the value Ch = (11/3) x 2.1621 = 7.927 
in Eqs. (3.41) and (3.68), which gives 

7 = 0.5562 - 199.2 Ca + O(c^) . (4.4) 

For the scalar box case discussed in Sec. (3.4) one has instead = 2.1621 and in this case one 
finds 

7 = 0.5562 - 54.8 Ca + 0{cl) . (4.5) 

As a comparison, we have also computed the exponent 7 for the case c/^ = in Eqs. (3.41) and 
(3.68). This corresponds to a case where the 0{h) correction to Spvac is entirely neglected, and 
one obtains 7 = 0.5562 — 0.703 + O(c^). Finally for the Newtonian (non-relativistic) treatment, 
described in Appendix A, one finds the much smaller correction 

7 = 0.5562 - 0.0142 Ca + 0{cl) . (4.6) 

Among these last expressions, the tensor box case is supposed to give ultimately the correct answer; 
the scalar box case only serves as a qualitative comparison, and the Ch = case is done to estimate 
independently the size of the correction coming from the ubiquitous 0{h) or terms (see for 

example the differential equation for the density perturbations S{t) in Eq. (3.75)). Note that the 
Ch = 0, scalar and tensor box results can be summarized into the slightly more general formula 

7 = 0.5562 - (0.703 + 25.04 Ch)ca + O(c^) . (4.7) 

showing again the overall importance of the Ch contribution to dp^ac in Eq. (3.41). This last term is 
responsible for the feedback of the metric fluctuations h on the vacuum density Spyaf, and pressure 
Spvac fluctuations. 

It should be emphasized here once again that all of the above results have been obtained by 
solving the differential equation for S{a), Eq. (3.92), with G{a) given in Eq. (3.11), and exponent 
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1p = 3/2i^ ~ 9/2 relevant for a matter dominated background universe. It is this last choice that 
needs to be critically analyzed, as it might give rise to a definite bias. Our value for so far 
reflects our choice of a matter dominated background. More general choices, such as an "effective" 
= 3(1 + w)/2u with and "effective" w, or even the use of Eq. (3.85), are in principle possible. 
Then, although Eq. (3.75) for 5{t) remains unchanged, Eq. (3.92) for 5{a) would have to be solved 
with new parameters. In the next section we will discuss a number of options which should allow one 
to increase on the accuracy of the above result, and in particular correct the possible shortcomings 
coming so far from the specific choice of the exponent 

4.1 Possible Physical Interpretation of the Results 

Looking at these last results (see also Figs. 1-4), they seem to indicate that (a) the correction due 
to the h (or l/2i/) terms in Eq. (3.41) and in the differential equation, Eq. (3.75), for 5{a) is rather 
large, and that (b) it is more than twice as large in the tensor box case than it is in the scalar box 
case. Furthermore they seem to suggest that (c) the Newtonian (non-relativistic) result, which does 
not contain a p^ac contribution, substantially underestimates the size of the quantum correction. 
To quantitatively estimate the actual size of the correction in the above expressions for the growth 
index parameter 7, and make some preliminary comparison to astrophysical observations, some 
additional information is needed. 

The first item is the coefficient cq ~ 33.3 in Eq. (2.28) as obtained from lattice gravity cal- 
culations of invariant correlation functions at fixed geodesic distance [19]. We have re-analyzed 
the results of [19] which involve rather large uncertainties for this particular quantity, nevertheless 
it would seem difficult to accommodate values for cq that are more than an order of magnitude 
smaller than the quoted value. A renewed more accurate lattice calculation of cq, obtained from 
the computation of invariant curvature correlation functions at fixed geodesic distance, would seem 
rather desirable at this point. 

The next item that is needed here is a quantitative estimate for the magnitude of the coefficient 
C(i in Eq. (3.11) in terms of q in Eq. (2.18), and therefore in terms of cq in the original Eq. (2.28). 
First of all one has ~ q, if is identified with a scale factor corresponding to a universe of size 
^; to a good approximation this corresponds to the universe "today" , with the relative scale factor 
customarily normalized at that time to a/ao = 1, although some large conversion factor might be 
hidden in this perhaps naive identification (see below). 
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Regarding the numerical value of the coefficient q itself, it was found in [13] that in Eq. (2.18) 
ct ~ 0.785 Co in the scalar box case, and a — 0.450 cq in the tensor box case. In both cases these 
estimates refer to values obtained from the zeroth order covariant effective field equations. In the 
following we will take for concreteness the average of the two values, thus ct ~ 0.618 cq. Then for 
all three covariant calculations recorded above Ca ~ 0.618 x 33.3 ^ 20.6, a rather large coefficient. 

From all of these considerations one would tend to get estimates for the growth parameter 
7 with rather large corrections! For example, in the tensor box case the corrections would add 
up to — 199. Ca = —199. X 0.618 x 33.3 = —4095.. Even in the Newtonian (non-relativistic) case, 
where the correction is found to be the smallest, the corresponding result appears to be quite large. 
In this last 2.7 Co (see Appendix A), so the correction to the index 7 becomes 

-0.0142 x 2.7 X 33.3 = -1.28. 

It would seem though that one should account somewhere for the fact that the largest galaxy 
clusters and superclustcrs studied today up to redshifts 2; ~ 1 extend for only about, at the very 
most, 1/20 the overall size of the visible universe. This would suggest then that the corresponding 
scale for the running coupling G{t) or G{a) in Eqs. (2.18) and (3.11) respectively, should be reduced 
by a suitable ratio of the two relevant length scales, one for the largest observed galaxy clusters or 
superclusters, and the second for the very large, cosmological scale ^ ~ \/^/\J^ ~ 1.51 x lO^^cm 
entering the expression for 5G{U) in Eqs. (2.3) and (2.28). This would dramatically reduce the 
magnitude of the quantum correction by as much as a factor of the order of (1/20)''''' = (1/20)^'^ 
1.398 X 10~^. When this correction factor is roughly taken into account, one obtains the more 
reasonable (and perhaps observationally more compatible) estimates for the tensor box case 

7 = 0.5562 - 0.0057 c^ + O(c^) . (4.8) 

and for the scalar box case 

7 = 0.5562 - 0.0016 c^ + O(c^) . (4.9) 

while in the non-relativistic (Newtonian) case one finds 7 0.5562 - 4.08 x 10"^ . In the tensor box 
case this would then amount to a slightly reduced value for the growth index 7 at these scales as 
compared to the constant G case, by as much as a few percent, which could perhaps be observable 
in the not too distant future. Of course, on larger scales the effects would be more significant, and 
somewhat bigger for larger values of il. 

A second possibility we will pursue here briefly is to consider a shortcoming, mentioned pre- 
viously, in the use of a{t) ~ ao(t/to)'^^^ in relating G(a) in Eq. (3.11) to G{t) in Eq. (2.18). In 
general, if w is not small, one should use instead Eq. (3.85) to relate the variable t to a{t). The 
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problem here is that, loosely speaking, for w ^ at least two itj's are involved, w = (matter) 
and w = —1 (A term). Unfortunately, this issue complicates considerably the problem of relating 
6G{t) to SG{a), and therefore the solution to the resulting differential equation for S(a). As a 
tractable approximation though, one should set instead a{t) - ao{t/to)W+^)^ and then use an 
"effective" value of w ^ —7/9, which would seem more appropriate for the final target value of 
CI Ri 0.25. For this choice one then obtains a significantly reduced power in Eq. (3.11), namely 
'ji, = 3(1 + w)/2u = 1. Furthermore, the resulting differential equation for (5(a), Eq. (3.92), is still 
relatively easy to solve, by the same methods used in the previous section. One now finds 

7 = 0.5562 - (0.92 + 7.70 Ch) Ca + O(c^) . (4.10) 

which should be compared to the previous result of Eq. (4.7). In particular for the tensor box case 
one has again Ch = 7.927, which can the be used to compare to the previous result of Eq. (4.4). 
Thus by reducing the value of ji, by about a factor of four, the Ca coefficient in the above expression 
has been reduced by about a factor of three, a significant change. 

After using this improved value for the power ^y^,, the problem of correcting for relative scales 
needs to be addressed again, in light of the corrected estimate for the growth exponent parameter 
of Eq. (4.10). Given this new choice for 7,^ = 1, on can now consider, for example, the types of 
galaxy clusters studied recently in [23, 24, 25], which typically involve comoving radii of '~ 8.5Mpc 
and viral radii of ~ lAMpc. For these one would obtain an approximate overall scale reduction 
factor of (1.4/4890)"^ ~ 2.9 x 10"'^. Note that in these units (Mpcs) the reference scale appearing 
in G(n) is of the order of ^ ~ 4890Mj9c. This would give for the tensor box (ch = 7.927 ) 
correction to the growth index 7 in Eq. (4.10) the more reasonable order of magnitude estimate 
—62. X 20.6 X 2.9 x 10~^ —0.37, and for 7 itself the reduced value would end up at ^ 0.19. Clearly 
at this point these should only be considered as rough order of magnitude estimates. 

Nevertheless this last case is suggestive of a trend, quite independently of the specific value of 
Ch and therefore of the overall numerical coefficient of the correction in Eq. (4.10): namely that the 
correction to the growth index parameter will increase close to linearly (for 7,^ close to one, as we 
have argued) in the size of the cluster. Consequently one expects that the deviations will increase 
tenfold in going from a cluster size of IMpc to one of WMpc, and a hundredfold in going from 
IMpc to lOOMpc. 

Finally another possible, and ultimately much more conservative, approach would be to take 
- at least for the time being - with some caution the rather large value for cq obtained from 
nonperturbative lattice quantum gravity calculations. One could then use instead the observational 
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bounds on x-ray studies of large galactic clusters at distance scales of up to about 1.4 — 8.5Mpc 
[24], namely 7 = 0.50 it 0.08, to constrain the value of the constant Ca at that scale, giving for 
example from Eq. (4.10) the bound Qj < 8 x 10~^ in the case of tensor box, and the much less 
stringent bound Ca < 0{l) for the Newtonian (non-relativistic) case of Eq. (4.6). 
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Figure 1. Illustration of the matter density contrast 6 (a) as a function of the scale factor a{t), in the fully 
relativistic treatment (tensor box) and for a given matter fraction 57 — 0.25, obtained from the solution of the 
density contrast equation of Ec^. (3.76), with G{a) given in Eq. (3.11) with j^, — 9/2 and for Ca — 0.001. In 
the case of a running G(n), one generally observes a slightly faster growth rate for later times, as compared 
to the solution for the case of constant G and with the same choice of fi, described by Eq. (3.87). 
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Figure 2. Illustration of the growth index parameter 7 of Eq. (4.3) as a function of the matter density 
fraction 57, computed in the Newtonian (non-relativistic) theory with a running G(a) given in Eq. (3.11), 
and obtained by solving Eq. (A. 36), here with with 71, = 9/2 and Ca — 0.01. For the specific choice of 
matter fraction 57 = 0.25, suggested by KCDM models, one then obtains the estimates for the growth index 
parameter given in Eq. (4.6). 
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Figure 3. Illustration of the growth index parameter 7 of Eq. (4.3) as a function of the matter density 
fraction fl, computed in the fully relativistic (tensor box) theory with a running G{a) as given in Eq. (3.11), 
and obtained by solving Eq. (3.75) with with ji, — 9/2 and Cq = 0.0003. For the specific choice of matter 
fraction ft — 0.25 one then obtains the estimates given for the tensor box in Eq. (4.4). Not surprisingly the 
deviations from the standard result for 7 become more visible for larger values of fl. 
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Figure 4. Qualitative comparison of the growth index parameters 7 of Eq. (4.3) as a function of the matter 
density fraction fl, computed first in the relativistic (tensor box) theory with a running G{a) and Ca — 0.0003, 
then in the Newtonian (non-relativistic) treatment also with a running G{a) and Cq = 0.01, both with 
7^ = 9/2, and finally compared to the usual treatment with constant G. In both cases the deviations from 
the standard result for 7 are most visible for larger values of il, corresponding to a greater matter fraction. 
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4.2 Density Perturbations in the Conformal Newtonian Gauge with G(n) 



In this section we will outline briefly what other avenues can be pursued to determine quanti- 
tatively and systematically the cosmological effects of a running G{0). The perturbed RW metric 
is well suited for discussing matter perturbations, but occasionally one finds it more convenient to 
use a different metric parametrization, such as the one derived from the conformal Newtonian (cN) 
gauge line element (see for example [31, 32], and references therein) 

dr^ = a'^{t){{l + 2ilj)dt'^ - {l-2(f>)Sijdx'dx^} (4.11) 

with conformal Newtonian potentials 4>{x.,t) and (?!)(x, f). In the simplest framework, the two 
potentials ip and (p give rise separately to Newton's equation for a point particle, and Poisson's 
equation, respectively 

x = -VV' = A7r Ga^ dp . (4.12) 

In this gauge, and in the absence of a G{D), the unperturbed equations are 

in the absence of spatial curvature {k = 0). In the presence of a running G these again need to 
be modified, in accordance with Eqs. (2.17), (2.16) and (2.18). A cosmological constant can be 
conveniently included in the p and p, with p\ = X/SttG = —p\. In this gauge scalar perturbations 
are characterized by Fourier modes V'(q;i) and (jii(q, i), and the first order Einstein field equations 
in the absence of G{p) read [31] 

a \ a J 

k^(^^+'^ip^ = A7rGa^{p + p)9 

a ^ ^ \ a a'^ I 3 3 

A:2(</,-^) = 12t: G {p + p) a (4.14) 

where the perfect fluid energy-momentum tensor is given to linear order in the perturbations 
dp = p — p and Sp = p — phy 

T% = -ip + Sp) 

r° = {p + p)vi = -T\ 

T) = {p + Sp)d'j + ^'j S\ = (4.15) 
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and one has allowed for an anisotropic shear perturbation S*^ to the perfect fluid form T^-. The 
two quantities 6 and a are commonly defined by 

{p + p)e = i y 5T] {p + p)a = -{kikj - (4.16) 

with = T*- — 5'^r|/3 the traceless component of T!^-. For a perfect fluid 6 is the divergence of 
the fluid velocity, 9 = ik^vj, with = dx^ /dt the small velocity of the fluid. The field equations 
imply, by consistency, the covariant energy momentum conservation law 

S = -(i + w)ie-34>)-3- (^-w)5 

a \dp ) 

e = _«(i_3^)^_^^e + ^^fe25-fcV + A:V (4.17) 
a 1 + w I + w dp 

and relate the matter fields 6, a and 9 to the metric perturbations (p and ip. where d is the matter 

density contrast S = 6p/p, and w is the equation of state parameter w = p/p. In General Relativity 

(f) = ip as long as there is no anisotropic stress, but in extended theories of gravity, such as the one 

described here, the relation between cf) and tp can become scale dependent. 

In the presence of a (?(□) the above equations need to be re-derived and amended, starting from 

the covariant field equations of Eq. (2.3) in the cN gauge of Eq. (4.11), with zeroth order modified 

field equations as in Eqs. (2.16) and (2.17), using the expansion for G{D) given in Eq. (3.24), but 

now in terms of the new cN gauge potentials (p and ijj. One key question is then the nature of 

the vacuum-polarization induced anisotropic shear perturbation correction S*^- appearing in the 

covariant effective field equations analogous to Eqs. (4.14), but derived with a 0(0). In particular 

one would expect the quantum correction to the energy momentum tensor appearing on the r.h.s. 

of Eq. (2.3) to contribute new terms to the last of Eqs. (4.14), which could then account for a 

non-zero stress a, and thus for a small deviation from the classical GR result for a perfect fluid, 

(p = tp. Naively one would expect ip/(p = 1 -|- 0{5G/Gq). An explicit calculation with (?(□) [35] 

gives 

^ = 1+1- ^ J 'iWvac = 1 + 1 - — — ^ 4.18 

(p) \ 2v[\^w)) Go \ 2uJ Go 

for u; = and Wvac = \- It is often customary (see e.g. [31, 32, 33, 34]) to parametrize deviations 
from General Relativity in terms of a slip function E and of the growth rate parameter 7 introduced 
previously. These two quantities are defined by 

V^{4) + ^P) = 3^VtH^5 1=^^H^ (4.19) 

logii 

with 5 the density contrast and / the density contrast exponent. Occasionally the parameter 
rj = ip/(p — 1 is introduced as well. In classical General Relativity ip/(p = l, r] = 0,'S = l and then 
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the growth exponent 7 ~ 0.55 for Q ~ 0.25. The calculations presented in the previous sections 
have already suggested to some extent what changes to expect for the exponent 7, which then 
leaves the problem of determining the structure of the S correction. In addition, the Newtonian 
(non-relativistic) calculation of Appendix A has determined, from the form of the modified Poisson 
equation, one of the relevant equations, namely the one for the potential 0. We plan to discuss 
these interesting questions in a future publication [35]. 



5 Conclusions 



In this paper we have attempted to systematically analyze the effects on matter density per- 
turbations of a running appearing in the original effective, non-local covariant field equations 
of Eq. (2.3). The specific form of G{D) in Eq. (2.1) is inspired by the non-perturbative treatment 
of covariant path integral quantum gravity, and follows from the existence of a non-trivial fixed 
point in G of the rcnormalization group in four dimensions. The resulting effective field equations 
are manifestly covariant, and in principle besides the genuinely non-perturbative scale ^ there are 
no adjustable parameters, since the coefficients (cq) and scaling dimensions (v) entering G{0) are, 
again in principle, calculable by systematic field theory and lattice methods ([6] and references 
therein) . 

The present work can be viewed in broad terms as consisting of two parts. In the first part we 
have systematically developed the general formalism necessary to deal with small matter density 
fluctuations in the presence of a running gravitational coupling G{U). Most, if not all, of the 
results in the first part have been formulated in a way that assumes as little as possible about 
specific aspects related to how exactly G does run with scale. Indeed many of the equations we 
have obtained are not restricted to v = i, and are found to be valid for a wide range of powers v 
and coefficients cq appearing for example in the original expression for G(n) as given in Eq. (2.28). 
Furthermore, the zeroth order (in the fluctuations) results of [13], on which the present work builds 
up, do not rely on any specific value for these parameters either, since the expressions obtained 
there follow from general properties of the covariant d'Alembertian and its powers, as they appear 
in G{n). In particular the flow in the vicinity of the ultraviolet flxed point could in principle allow 
for Co being either negative (gravitational screening) or positive (gravitational anti-screening), and 
both cases could in principle be described by the results obtained above, for example for the growth 
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index / and the growth index parameter 7. It is only the latter option though that is favored by 
studies of non-perturbative Euclidean lattice gravity (the weak coupling phase is unstable and 
found to describe a collapsed degenerate two-dimensional spacetime), hence the choice here to 
discuss primarily this last case. But in principle the fact remains that the sign of cq will ultimately 
determine the direction of the corrections given above, which could eventually become constrained 
by observation. In the end the only result that is extensively used in the first part is the result 
of [13] that Wyac = 5) apart from the fact that we choose to restrict our attention from the very 
beginning primarily to the non-relativistic matter case w = 0, and to the large wavelength limit 
q ^ 0. Later on it was found that for sufficiently slowly varying backgrounds the result 
preserved also to first order in the perturbations, which seems to suggest some level of consistency 
in the treatment of the field perturbations. 

In spite of the non-locality of the original effective field equations in Eq. (2.3), one finds quite in 
general that small perturbations can be treated, in a first approximations, in terms of local terms, 
described by quantities p^ac and p^ac as they appear in the effective description of T^"*^ in terms of 
a perfect fluid. The latter should then be regarded as the leading term in a derivative expansion of 
the non-local contribution to the effective field equations, as they apply here to the rather specific 
case of the FRW background. Under the physically motivated assumption of a comparatively 
slowly varying (both in space and time) background, it is then possible to obtain a complete and 
consistent set of effective field equations, describing small perturbations for the metric trace and 
matter modes (Eqs. (3.69), (3.72), (3.73) and (3.74)). From these a single equation for the matter 
density contrast is eventually obtained, Eq. (3.75), which is the main result of this work. The only 
input needed in this last equation is dG{t), the zeroth order (in the fluctuations) running of G as 
written in Eq. (2.18), with given more or less known parameters u and Cj. The corresponding result 
in the Newtonian (non-relativistic) treatment is obtained in Appendix A, leading to Eq. (A.30). 

The next step was a translation of the equation for the density contrast S{t) into the correspond- 
ing equation for S{a), involving a related running coupling G{a), instead of the original G{t). Since 
in general the transformation from one variable to the other is not entirely trivial, some simplifica- 
tion had to be assumed, i.e. that the quantum correction in G{a) can be written as a power, with 
an exponent 7^/, a choice that could in the future be relaxed as part of a broader more systematic 
investigation. Subsequently a solution for the differential equation for 6{a) was obtained, leading to 
expressions for the growth index /(a) and for the growth index parameter 7. A number of general 
features can be observed, the first one being the fact that generally the correction to the growth 
index parameter 7 is found to be negative, indicating a less steep rise of / with ft. 
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The second part of the paper describes a number of attempts to provide a semi-quantitative 
estimate for the corrections obtained, in order to see whether these corrections could be related 
in some way to current astrophysical observations. In order to do so, one needs to adapt the 
theoretical calculation for the growth index parameter 7 to the kind of observational data available 
from the study of large galactic clusters. This requires, as expected, a careful consideration of 
the relative length scales that come into play. On the one hand, one length scale is given by the 
size of the largest clusters reached by observation, typically of the order of a few Mpcs. On the 
other hand it should involve the absolute reference scale given by ^ = ^/3/X 2± 4890Mpc. The 
comparison between theory and observation would then seem straightforward, were it not for the 
fact that this ratio generally comes in to a certain power, whose detailed knowledge is necessary 
in order to eventually reduce the quantitative uncertainties. Eventually these could be bracketed 
by a more systematic study of the solutions to the 6{a) equation, and the corresponding growth 
exponents 7. We are referring here in particular to a study of the sensitivity of the results to the 
specific choices of the exponent 7,^, appearing in SG{a) and determined in part by the relationship 
between the variables t and a{t), which we discussed earlier. In addition, there is still perhaps a 
certain level of uncertainty in the actual coefficients cq and q entering the theoretical predictions, 
which we have also described above in some detail. The latter could be reduced further by improved 
non-perturbative lattice computations. Nevertheless, the value of the present calculations lies in 
our opinion in the fact that so far a discernible trend seems to emerge from the results. The trend 
we have found seems to suggest that the correction to the growth exponent 7 is initially rather 
small for small clusters, negative in sign, and then slowly increasing in magnitude, close to linearly 
with scale. 

It is clear that the effects discussed in this paper are only relevant for very large scales, much 
bigger than those usually considered, and well constrained, by laboratory, solar or galactic dynamics 
tests [1, 20, 21, 22]. Furthermore the effects we have described here are quite different from what 
one would expect in f{R) theories [26, 27], which also tend to predict some level of deviation from 
classical GR in the growth exponents [28, 29, 30]. Future more accurate astrophysical observations 
might make it possible to see the difference in the predictions of various models [33, 34, 36, 37, 38]. 

Acknowledgements 

One of the authors (HWH) wishes to thank Thibault Damour and Gabriele Veneziano for in- 
spiration and discussions leading to the present work, and Alexey Vikhlinin for correspondence 
regarding astrophysical measurements of structure growth indices. He also wishes to thank Her- 



42 



mann Nicolai and the Max Planck Institut fiir Gravitationsphysik (Albert-Einstcin-Institut) in 
Potsdam for warm hospitality. Thanks also go to Salvatorc Capozziello and Paolo Serra for pro- 
viding us with a number of useful references. The work of HWH was supported in part by the Max 
Planck Gesellschaft zur Forderung der Wissenschaften, and by the University of California. The 
work of RT was supported in part by a DoE GAANN student fellowship. 



43 



Appendix 



A Non-Relativistic (Newtonian) Treatment of Matter Density Per- 
turbations 



In this section we discuss the Newtonian theory of small matter fluctuations, first by recalling 
the relevant equations in the usual treatment, and then by presenting what changes need to be 
implemented in order to account for the running of G. Later these equations will be solved, so that 

a comparison can be made with the results in the absence of a running G. 

When discussing a nonrelativistic Hubble flow it is customary to define coordinates in the 
following way 

r dr a .. ^, 

^=—rT v = — = -r (A.l) 

a{t) dt a ^ ' 

where x is attached to the comoving frame, while r is the flat Minkowski space coordinate, such 
that in the comoving frame x one has, by construction, dx/dt = 0. 

In the following some simplification will arise due to the fact that we shall consider a non- 
relativistic fluid with the negligible pressure, p ~ or w; = 0. The relevant equations are then 
the continuity equation, the Eulcr equation and the gravitational field equations. These will be 
listed below to zeroth and first order in the matter density (p), pressure (p), velocity field (v) and 
gravitational field g. 



A.l Newtonian Treatment Without the Running of G 

After decomposing the fields into a background and a fluctuation contribution, p = p + 5p^ 
p = p + 6p, and v = v + 6v, one obtains from the continuity equation, to zeroth and first order 

respectively, 

p + V-(^v) = 6p + 3-6p+- (r-V) 6p + pV -Sv = . (A.2) 

a a 

When the effect of the Hubble flow is included, i.e., Eq. (A.l), the above zeroth order equation 
reduces to 

m+^^^m =0 (A.3) 
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with solution p{t) = pQ (ao/a(t))^, where po and are the two integration constants corresponding 
to the present matter density and to the present scale factor (usually taken to be ao = 1). We note 
here that Eq. (A. 3), and hence Eq. (3.14), will continue to hold for a running G, as these equations 
are derived from the kinematics and the continuity equations in the RW background metric given 
in Eq. (A. 2), whose is not affected by the running of G — >■ G{n). 

To zeroth and first order in the fluctuations the Euler equations for a fluid in the RW background 
are given respectively by 

V + (v • V) V = g S'y + - S\- + - (r ■ V) 5v = ~V 6p + 6g . (A.4) 

a a p 

Finally the gravitational field equations are given to zeroth and first order in the fluctuations by 

Vxg = V-g = -47rGop (A.5) 

V X 5g = V • (5g = - 47r Go (5/9 (A.6) 

incorporating Gauss' law and the constraint that the gravitational fields are longitudinal. Only 

the last set of equations contain the gravitational constant G. Hence, in the framework of the 
Newtonian treatment, the modification of a running G G(0) only affects the gravitational 
Poisson equation. 

It is customary at this stage to introduce Fourier components of the fluctuations, and write 



Sp{r,t) = 6pq{t) exp 



z r • q 



(A.7) 



ait) 

and similarly for (5v, Sg , and possibly Sp. For an adiabatic fluctuation one can also set 6p = Vg Sp, 
with Vg the speed of sound. 

Then to first order in the fluctuations the continuity equation, Euler equation and the gravita- 
tional field equations take on the form, for each mode q, 

^''^^'^ + ^ W) ^'""^'^ ^ '""'aitf^ '^'^ = ° ^^-^^ 
^v,(t) + ||^v,(0 = -^|^^P.(t) + ^g,(t) (A.9) 

5gq(i) = ^ a{t) Go dp^it) . (A.IO) 

Subsequent elimination of the gravitational and velocity fields then leads to a single second order 
differential equation for the matter density contrast Sq{t) = Spq{t)/p{t) describing the physics of 
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compressional modes: 

S^{t) + 2 ^ + - An Go S^{t) = . (A.ll) 

In the limit of very long wavelength fluctuations, q — > 0, the above equation simplifies to 

6{t) + 2 ^ 5{t) - Att Go p{t) 5{t) = . (A.12) 
a[t) 

A solution can then be found, using p{t) = l/QirGt^ and a{t)/a{t) = H{t) = 2/3t, such that 
the general form for 6{t) is given by a linear combination of cither ~ t"^^'^ or ~ t^^. The latter 
corresponds to a decaying (as opposed to growing) solution and is usually discarded, giving finally 
the standard Newtonian result 6{a) oc a. Wc note here that the above non-relativistic equation and 
solution applies to the case of non-relativistic matter only; in particular it excludes the presence of 
a cosmological constant. 



A. 2 Newtonian Treatment with Running G(n) 



The next step is a modification of the non-relativistic equations in Eqs. (A. 2), (A. 4), (A. 5) and 
(A. 6) to incorporate a suitable running of G. Since only the latter set of equations, Eqs. (A. 5) 
and (A. 6), contain G it is only these that need to be suitably modified. In the presence of a 
scale-dependent coupling one has 

(5g = -V(5(/) (A.13) 
with the perturbing potential 64> given by a solution to Poisson's equation 



V^S(f){r, t) = -V- (5g(r, t) = 47rG(n) 6p{r, t) 



(A.14) 



and (?(□) given in Eq. (2.28). Following Eq. (A. 7), as it applies here to 6g and dp, we Fourier trans- 
form the spatial components of the above Poisson equation, which requires the Fourier transform 
of G{n) as obtained from Eq. (2.28), namely 



As a result the gravitational field perturbation is of the form 



+ ... 



4:7r i 

<^gq(i) = — o-oW • exp 



zr • q 



ait) 



G(q2,9|) (5pq(i)exp 



zr • q 



a{t) 



(A.15) 



(A.16) 
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Since we are mainly interested in the long wavelength limit, it suffices here to evaluate the above 
expression in the limit q — >^ 0, 



q 

47riq 



1-^ + ... 



ait) 



a{t) 

G{ci-\d^)5p^{t)- 



a{t) 



■G{e,d^)6p^{t) + G{ci\d^)Sp^{t) 



+ ... 



(A.17) 

and for q = only the first term survives. Furthermore, when (?(□) = G(q^, df) acts on a function 
of t which we will assume here is of the form of a power {e.g., t", with the power a a number of 
order one) one obtains 

G(q2, 5|) • r ^ G(t) • i° . (A. 18) 

Here the running coupling G{t) is given by the expression in Eq. (2.18), with to = ^, and the 
coefficient 

r(i + a) 



r(l + a + l/i/j 

Thus for example for a = —4/3 (the standard Newtonian result for matter density perturbations) 
one has q = (27/10) cq; in the following it will be safe to assume that the coefficient ct in Eq. (2.18) 
is a number of the same order of magnitude as the original cq in Eq. (2.28). 

Consequently, when acting on a density perturbation (5pq(t) in the form of a power law in t, to 
leading order in q one obtains simply 



5gq(0 = ^«(i)G0 



l + Ct 



+ ... 



5p^{t) . 



(A.20) 



This last result can be compared with Eq. (A. 10) for the case of a constant G. 

As stated previously, the continuity equation for the fluctuations, Eq. (A. 8), and the corre- 
sponding Eulcr equation for the fluctuations, Eq. (A. 9), arc not modified by the presence of a 
running G(n), as given in Eqs. (A. 16) and (A.20). To solve the resulting equations of motion 
for the fluctuations, it is now customary to decompose the velocity perturbation 5^r into parts 
perpendicular and parallel to q 



(^Vq(t) = (5vq_L(i) + iqeq(t) 



with 



q - (5vq_ 







i q • (5vq 



Cq r, 



The fractional change in the matter density 5 is then defined as 



Pit) 



(A.21) 

(A.22) 
(A.23) 
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With the above decomposition of the velocity field 6v and the expression for the density contrast 
5 inserted into the first order continuity equation, Eq. (A. 8), one obtains the unmodified result 

^2 



ait) 



(A.24) 



so that there is no change in the relationship between S and e when G G{D). In turn the Euler 
equation for the fluctuation, Eq. (A.9), now becomes the two sets of equations 



Re : 6vq±{t) + -dvq±{t) = 



(A.25) 



with the gravitational field fluctuation Sgq now given by the expression in Eq. (A. 16). Prom the 
real part (corresponding to rotational modes) one concludes 



(A.26) 



which is of the same form as in the case of a constant G. From the imaginary part (corresponding 
to compressional modes) in Eq. (A.25) one obtains, using Eq. (A.24), 



2 

<5q(t) + 2 ^ jq(t) + ^vl <^q(t) " CXp 



zr ■ q 



a{t) 



G(q2,a|) ( exp 



ir ■ q 



p(t)<5q(t) =0. (A.27) 



The latter can be recast into the slightly simpler form 



a i 



6^{t) + 2 - 6^{t) +[^v^^-47rg (q^ d't) 1 6^{t) = 



(A.28) 



by defining a modified source term 



a(q^5^)^ 



1 



exp 



— ^r • q 



a{t) 



G(q^a/) exp 



ir ■ q 



ait) 



ms^it))] 



In the limit q — > one obtains immediately 



5it) + 2 ^ 6it) - 47r Git) pit) 5(t) = . 



(A.29) 



(A.30) 



The last two equations can now be compared with the corresponding results for a constant G, given 
in Eqs. (A. 11) and (A.12). 



A. 3 Computation of the Non-Relativistic (Newtonian) Growth Index with G(n) 
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The next step requires a solution of the differential equation for the density perturbations (Jq(t), 
in the Newtonian approximation and in the limit q — t- 0, as in Eq. (A. 30). It is convenient and 
customary at this point to change variables from t to the scale factor a{t), so that 6q{t) — >■ 5q{a) = 
Sq ■ 5{a). Prom Eq. (3.90) one has 

d5{a) 



m 
m 



aH{a) 
a^H\a) 



da 
dlnHja) 
da 



+ 



(A.31) 



Here H{a) is defined as the Hubble "constant" H{a) = a{t)/a{t), as it appears in the equations of 
motion for a background FLRW geometry 



H{a) = JfG{a)p{a) + ^, 



but with a running Newton's constant G{a) (see Eq. (2.18)) 

SG{a 



G{a) = Go 



1 + 



1 + Ca 



ao 



+ ... 



(A.32) 



(A.33) 



Here the index is ji, = S/2u, since from Eq. (2.18) one has for non relativistic matter a{t)/ao ^ 
{t/tof/^. In the above expression Cq ~ Q if ao is identified with a scale factor corresponding 
to a universe of size ^; to a good approximation this corresponds to the universe "today", with 
the relative scale factor customarily normalized to a/ao = 1. As a consequence, the constant Ca 
in Eq. (A.33) can be taken to be of the same order as the constant co appearing in the original 
expressions for (?(□) in Eqs. (2.5) and (2.28). Note also that, by the use of Eq. (A.32) for the 
scale factor, wc have allowed for a non- vanishing cosmological constant in our otherwise Newtonian 
(non-relativistic) treatment. 

After these substitutions one finally obtains the differential equation for the matter density 
contrast, Eq. (A.30), in the variable a{t) 



(f5{a) 



+ 



d\iiH{a) 3 
da a 



d5{a) 47rG'(«) p{a) 
da 



S{a) = . 



(A.34) 



a2 H^{a) 

Note that in order to compute the leading, in 5G{a)/Go, correction to the density contrast S{a), 
one only needs p{a) to lowest order as given in Eq. (3.14), and H{a) as given in Eq. (A.32). 
With the aid of the parameter 6 (see Eq. (3.86)) 

1-n 



e = 



(A.35) 
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where ft is the matter density fraction and 1 — 17 the cosmological constant fraction as measured 
today, one obtains the following differential equation for the density contrast 5{a) 

3(l+2a3g) / lua"' + {hu-l) a^-g \ d5_ _ 3 / a^- g \ 

5a2 + 2a(l+a3^) r^""" {l+a^ 6) {\+2a^ 9) j da 2a^ {l+a^ 0)\-^'''' l+a^ 6 ] 

(A.36) 

for a reference scale ao = 1; the latter can always be re-introduced later by the trivial replacement 
a — >■ a/ao- 

Without a scale-dependent G (cq = in Eq. (A.33)), the growing solution to the above equation 
is given by 

5o{a) oc a • 2^1 Q, 1; y ; -a^ ^) (A.37) 

where 2^1 is the Gauss hypergeometric function. To evaluate the correction to 5o(a) coming from 
the terms proportional to Cq one sets 

(5(a) oc a-2i^i (^,l;^;-aH^ [l + c„^(a)] , (A.38) 

then inserts the resulting expression in Eq. (A.36), and finally expands the resulting expression 
to lowest order in to find the correction J-'(a). The resulting differential equation can then be 
solved for J- (a), giving the density contrast 6{a) as a function of the two parameters (ji, and or 
9 = (1 — i7)/il) appearing in Eq. (A.36). In the following we will focus on the specific choice = ^ 
obtained from the lattice theory of gravity [8], which leads to the G{a) exponent 7,^ = ^ = 9/2. It 
is customary at this point to define the growth index /(a) = ^g*^,^^"^ and the related growth index 
parameter 7 via 7 = . Then the solution to Eq. (A.36) gives an explicit expression for the 

a=ao 

growth index 7 parameter, as a function of the matter fraction il. 

Based on observational constraints, one is mostly interested in the case Q ^ 0.25, therefore 
in the following we will limit our discussion to this choice only. In the absence of a running G 
(G Go, thus Ca = 0) one has f{a = ao) = 0.4625 and 7 = 0.5562 for Q = 0.25 [17]. On the other 
hand when the running of G is taken into account one finds from the solution to Eq. (A.36) for the 
growth index parameter 7 at = 0.25 

7 = 0.5562 - 0.0142 Ca + O(c^) . (A.39) 

In the end it would seem therefore that at least in the Newtonian treatment the correction comes 
out rather small. Note that both the Newtonian and the relativistic treatment, described in the 
main body, give a negative sign for the correction arising from the running of G. 
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To estimate quantitatively the actual size of the correction in Eq. (A. 39) one needs an estimate 
for the coefficient cq ~ 33.3 in Eq. (2.28), as obtained from the lattice gravity calculations of 
invariant correlation functions at fixed geodesic distance [19]. In addition one uses the fact that 
Ca ct ^ 2.7 Co (see the previous discussion related to Eq. (A. 33)). From this one would then get 
the estimate 7 = 0.5562 — 1.28 on the largest scales, which looks like a significant 0(1) correction 
to 7. 
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